Determination of the angular momentum 
distribution of supernovae from gravitational wave 

observations 

'66 

(3 '. K Hayamai, S Desai^, K Kotake^, S D Mohantyi, 

CSJ ' M Rakhmanov^, T Summer scales'', S Yoshida^ 

^n. ^Thc University of Texas at Brownsville, Brownsville, TX, 78520, US 

r^ ' ^ Center for Gravitational Wave Physics, Pennsylvania State University, 

^ • University Park, PA 16802, US 

^v( ' ■^Division of Theoretical Astronomy/Center for Computational Astrophysics, 

National Astronomical Observatory of Japan, 2-21-1, Osawa, Mitaka, Tokyo, 
181-8588, Japan 

■^Department of Physics, Andrews University, Berrien Springs, MI 49104, US 
* Southeastern Louisiana University, Hammond, LA 70402, US 



cn 









in 



E-mail: kazuhlro.hayamaaiigo.org 



Abstract. Significant progress has been made in the development of an 
international network of gravitational wave detectors, such as TAMA300, LIGO, 
VIRGO, and GEO600. For these detectors, one of the most promising sources of 
fv^ , gravitational waves are core collapse supernovae especially in our Galaxy. Recent 

►-^ , simulations of core collapse supernovae, rigorously carried out by various groups, 

show that the features of the waveforms are determined by the rotational profiles 
of the core, such as the rotation rate and the degree of the differential rotation 
prior to core-collapse. Specifically, it has been predicted that the sign of the 
second largest peak in the gravitational wave strain signal is negative if the core 
rotates cylindrically with strong differential rotation. The sign of the second peak 
^^, ■ could be a nice indicator that provides us with information about the angular 

f^^ ' momentum distribution of the core, unseen without gravitational wave signals. 

ryr\ ' Here we present a data analysis procedure aiming at the detection of the second 

^»^ ' peak using a coherent network analysis and estimate the detection efficiency when 

a supernova is at the sky location of the Galactic center. The simulations showed 
we were able to determine the sign of the second peak under an idealized condition 
. ■ of a network of gravitational wave detectors if a supernova occurs at the Galactic 

J^ ' center. 

PACS numbers: 04.80.Nn, 07.05.Kf, 95.55.Ym, 95.85.Nv 

1. Introduction 

Current estimates of Galactic core-collapse supernovae rates predict one core-collapse 
supernova event every ~ 40 ± 10 year [TJ [2]. When a massive star (>^ IOMq) 
[3] undergoes a core-collapse supernova explosion in our Galactic center, copious 
numbers of neutrinos are produced, some of which may be detected on the earth. 
Such supernova neutrinos will carry valuable information from deep inside the core. In 
fact, the detection of neutrinos from SN1987A (albeit in the Large Magellanic Cloud) 
paved the way for Neutrino Astronomy, an alternative to conventional astronomy by 
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electromagnetic waves [JJ [SJ [S] . Core-collapse supernovae are now poised to give birth 
to yet another astronomy, Gravitational-Wave Astronomy. Currently long-baseline 
laser interferometers such as LIGO, VIRGO, GEO600, and TAMA300 are operational 
[71 [HI [HI [ini- For these detectors, core-collapse supernovae have been proposed as 
one of the most plausible sources of gravitational waves (GWs) (see, e.g., [Til HI] for 
reviews). 

If gravitational collapse of the supernova core proceeds spherically, no GWs can be 
emitted. To produce GWs, the gravitational collapse should proceed aspherically and 
dynamically. Observational evidence gathered over the last few decades has pointed 
towards core-collapse supernovae indeed being generally aspherieal [131 [HI [H] ■ The 
most unequivocal example is SN1987A. Recent HST images of SN1987A are directly 
showing that the expanding envelope is elliptical with the long axis aligned with 
the rotation axis inferred from the ring. The aspect ratio and position angle of the 
symmetry axis are consistent with those predicted earlier from the observations of 
speckle and linear polarization. What is more, the linear polarization became greater 
as time passed [Ml [16] , a fact which has been used to argue that the central engine of 
the explosion is responsible for the non- sphericity [T71 [18] . From a theoretical point of 
view, clarifying what makes the dynamics of the core deviate from spherical symmetry 
is essential in understanding the GW emission mechanism. 

It is well known that stars are generally rotating [TH] (see [501 [H] for recent 
controversy about the angular momentum transfer problem during stellar evolutions). 
This stellar rotation has been long supposed to play an important role in the 
GW emission from core-collapse supernovae, because the large-scale asphericities at 
core bounce induced by rotation can convert part of the gravitational energy into 
the form of GWs. Such GW signals are often simply called bounce signals (see 
[22l [23l [24l [25l [26] for other sources of GWs from the postbounce phase such as 
convection and anisotropic neutrino radiations). Just recently, several investigations 
have been carried out that explore how astrophysical information can be obtained from 
supernovae GW waveforms [271 [28]. Hayama (2004) [27] studied how accurately the 
waveforms associated with various supernova models may be estimated. Summerscales 
et al.(2008) [55] developed a maximum entropy approaches to extract the waveforms 
from noisy observations and discover what source information they contain. Such 
data analysis methods are not only indispensable for the first direct GW observation 
to confirm Einstein's theory of General Relativity, but also for clarifying supernova 
physics itself; such as the long-veiled explosion mechanism, hidden deep inside the 
core. 

Theoretically, much attention has been paid to the bounce signals in the context 
of rotational (e.g., [29 1 130 1 [31 ] [32 1 133 1 [34] and references therein) and magneotrotational 
core-collapse [35l [36l [37] with different levels of sophistication employed in the 
treatment of the equations of state (EOSs), neutrino transport, and general relativity 
(see [T5] for a review). Waveforms obtained in numerical simulations are often 
categorized to the three types, namely types I, II, and III. Type I waveforms are 
distinguished by a large amplitude peak at core bounce and subsequent damping ring- 
down oscillations. Type I waveforms are obtained when the initial angular momentum 
is small, which leads to core bounce near nuclear densities. As the angular momentum 
of the core becomes larger with stronger differential rotation, the waveforms transit 
to type II waveforms, which have several distinct peaks caused by multiple bounces 
supported by centrifugal forces. Type III, which are obtained for an extremely soft 
equation of state, are now considered to be unrealistic. 
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Figure 1. The figures represent waveforms for various core rotational models 1311 . 
as discussed in the text. Top left, top right and bottom left plots in the figure 
show the waveforms of the model MCS, RCS, SCS supernova signals with the 
optimal inclination. The bottom right plot is a signal from a model with a shell 
type rotation law and long differential rotation which moans a small amount of 
differential rotation. 



Through inspection of some of the elaborate numerical simulations cited above, 
it was suggested that the sign of the second largest peak becomes generically negative 
when the core rotates cylindrically with strong differential rotation [3T]. Looking 
carefully, this feature also can be seen in the type II waveforms obtained by other 
simulations ([3H [351 ISZ])- If we can observe the sign of the second peak, we could 
obtain information about the angular momentum distribution of an evolved massive 
star. The absolute amplitude of the second peak is within the sensitivity of initial 
LIGO for such a source at the galactic center. By using world-wide detector networks, 
the gravitational waveforms h^, hx can be reconstructed. In this paper, we look at 
how we might detect the second peak by using current LIGO and VIRGO. We propose 
a method to recover gravitational waveforms from supernova core collapse by using the 
RIDGE coherent network analysis [39|, which takes full advantage of multiple detectors. 
By carrying out Monte Carlo simulations, we estimate the detection efficiency of 
determining the correct sign of the second peak using an idealized LIGO- VIRGO 
network. 



2. Second peak of gravitational v^raves from core-collapse supernovae 

We take the waveforms from Kotake et al. [nU , in which a series of two-dimensional 
Newtonian hydrodynamic simulations of the rotational collapse of a supernova core 
was done using a realistic equation of state and a treatment of neutrino transport via 
the leakage scheme. Since we know little about the angular momentum distributions 
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in the core of evolved massive stars, their simulations assumed two possible rotation 
laws, a shell-type and a cylindrical rotation law. Twelve models of gravitational waves 
were computed by changing the combination of total angular momentum: S (slow), 
M (moderate), R (rapid), the rotation law: S (shell-type), C (cylindrical), and the 
degree of differential rotation: L (long), S (short), which indicates the distance from 
the rotational axis for the cylindrical rotation or the distance from the radius for the 
shell-type rotation. From the simulations, they found that the sign of the peak with 
the second largest absolute amplitude is negative for models with strong differential 
rotation and with a cylindrical rotational profile in the core (model MCS, RCS, SCS) 
and positive for the others. Figure [T] shows the waveforms for model MCS (top left 
plot), RCS (top right plot), SCS (bottom left plot), and MSL (bottom right plot). 
This observation implies that we can obtain information about the angular momentum 
distribution of massive evolved stars if we can detect the signs of the second peaks. 
This feature of the second peak also appears in simulations by other groups working on 
simulations of core-collapse supernovae (See figure 6 in [3Dj, figure 11 in [32], figure 10 
in [38], figure 8 in [37]. [|) 

We briefly describe how the sign of the peaks is determined. Generally, the sign 
of a gravitational waveform becomes negative right before core bounce because the 
oblateness of the core increases rapidly. After the bounce, the sign becomes positive 
as the oblateness decreases due to the shock propagating from the inner core surface 
outwards along the rotational axis. This explains the first largest peak amplitude and 
origin of the negative sign of the first peak. After the first peak, the core oscillates 
with decaying absolute amplitude because of its inertia, known as a ring-down phase. 
In the ring-down phase, the absolute amplitudes of the subsequent peaks are much 
smaller than the first peak described above. After the first peak, the core expands and 
then begins to shrink again. The core's change from expansion to contraction causes 
the positive peak in the waveform, which we call peak 1. Afterwards the amplitude 
begins to decrease until the core rebounds again. The rebound causes the subsequent 
negative peak in the waveform, which we call peak 2. The sign of the second largest 
peak is determined by which absolute amplitude [peak 1 or peak 2) is larger. This 
is sensitive to the degree of differential rotation of the core. For stronger differential 
rotation, the centrifugal force becomes so strong that the gravitational contraction of 
the core becomes weak. This leads to a weaker core-bounce than for smaller differential 
rotation, given the same initial rotation rate. Thus peak 1 becomes smaller as the 
differential rotation becomes stronger. Moreover, the influence of differential rotation 
can become more prominent for models with an initial cylindrical rotation profile, 
which is expected to be a natural angular momentum distribution inside the core [19] . 
These factors create a tendency for the absolute amplitude of the negative peak [peak 
2) to be larger than that of the positive peak {peak 1) for stronger differential rotation 
with the cylindrical rotation profile. 

3. RIDGE coherent netw^ork analysis 

RIDGE is a coherent network analysis pipeline described in [39] . The pipeline consists 
of two main components, which are data conditioning and generation of detection 
statistics. 

X However, according to 1401 . the type II waveforms tend to transform to the type I when the general 
relativity and deloptonization effects are taken into account. 
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Figure 2. The detector antenna function of HI, LI, VI for the Galactic center 
around 3:00 am, December 25th in 2007. Here we take as the Galactic center, the 
sky location at RA = 16.3 hr, DEC = —15.8°. The root mean squared detection 
response of /i_|_ is 0.57. 



The aim of the data conditioning is to whiten the data, so as to remove frequency 
dependence and any instrumental artifacts in the data. In RIDGE, the data are 
whitened by estimating the noise floor using a running median j41j . The whitening step 
above is followed by a line estimation and removal step by using a method described in 
[32] ■ The resulting conditioned data are passed on to the next step consisting of the 
generation of a detection statistic. The basic algorithm implemented for the coherent 
network analysis in RIDGE is Tikhonov regularized maximum likelihood, described in 
[43] . Recent studies O |43l [45] show that the inverse problem of a response matrix 
of a detector network becomes an ill-posed one and the resulting variance of the 
solution becomes large. This comes from the rank deficiency of the detector response 
matrix. The strength of the rank deficiency depends on the sky location, and therefore, 
plus or cross-polarized gravitational wave signals from some directions on the sky 
become too noisy. In RIDGE, we reduce this ill-posed problem by posing the Tikhonov 
regulator, which is a function of the sky location. The input to the algorithm is a 
set of equal length, conditioned data segments from the detectors in a given network. 
The output, for a given sky location 6 and (j), is the value of the likelihood of the data 
maximized over all possible /i+ and hx waveforms with durations less than or equal 
to the data segments. The maximum likelihood values are obtained as a function of 
and 4> ~ this two dimensional output, S{6, (jj), is called a skymap. We can construct a 
detection statistic called the radial distance which uses the entire skymap and scales 
the maximum likelihood by the mean location of the same quantities in the absence 
of a signal. A detailed description of the detection statistic is given in [39] . 

4. Simulations and results 



We carried out Monte Carlo simulations to study the detectability of the second peak 
of a supernova signal using RIDGE. Gravitational waves used here are "-I-" mode 
polarization waves and the models used were the MCS, RCS, SCS from [3T]. We 
assume the rotational axis of the core is perpendicular to the direction to the Earth 
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Figure 3. Left figure shows the peaks and their regions. Right plot shows the 
fiow chart for the detection of the second peak. First, wo find the second peak and 
the third peak in a gravitational wave. When the second peak of the reconstructed 
waveform is within 5 ms of the true second peak, the algorithm checks whether 
the third peak is within 5 ms of the true third peak. If this condition is satisfied, 
we declare that the detection of the sign of the second peak is successful. 



from the supernova so that gravitational waves are most strongly emitted towards 
Earth. We assume the rotational axis of the core is perpendicular to the direction to 
the earth from the supernova so that the gravitational waves are most strongly emitted. 
The location of the source was assumed to be that of the galactic center. The network 
consisted of the three LIGO detectors (Hanford/Livingston) and the VIRGO detector. 
All data streams are white Gaussian noise with the same sensitivities so that we focus 
on the performance of the signal reconstruction of the coherent network analysis. 
Therefore we need to add data conditioning for realistic analysis. We assumed that 
the detection of the supernova signal occurred 3:00 am, December 25th in 2007, i.e. 
December 25th 2007 03:00. This particular event time was chosen so both the VIRGO 
and LIGO detectors have reasonable sensitivities to this signal. Figure [2] shows the 
detector response for the supernova signal. The X-axis represents time and the Y-axis 
shows the sum-squared detector response. Simulated data streams are generated at a 
sampling frequency of 16384 Hz. First, the pipeline down-sampled all the data streams 
by applying anti-aliasing filters, resulting in a sampling frequency of 2048 Hz. In regard 
to the time-frequency granularity, the single tiles are {Sf,St) = (1024 Hz, 0.49 ms). 
The supernova signal was injected at every other second, and, in all, 200 signals 
were injected. The data was passed on to the RIDGE algorithm. In RIDGE, the data 
conditioning step was omitted since the noise was white Gaussian. The data was 
passed on to the step of the regularized coherent network analysis, and then, the 
estimated polarization waveforms /i+ were reconstructed. Picking up the reconstructed 
data around the injected signals we find the second largest peak of absolute amplitude. 
If the positions of the second peak and the third peak are within 5 ms of the center 
of the true peak positions and the sign is negative, we argue that the detection of 
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Figure 4. Detection probabilities of the second peak of the model MCS supernova 
at the galactic center is the left upper plot, the model RCS is the right upper plot, 
and the model SCS is the bottom plot in the figure. The detection probability 
is defined as a ratio of the number of the "-" sign detection divided by the total 
number of trials which is equal to 200. The signal-to-noise ratio on the X-axis is 
averaged over all detectors. 



the sign of the second peak suecceds, otherwise the detection fails. This efficiency is 
characterized by a detection probabihty which is defined as the ratio of the number of 
the detections to the total number of the injected signals. This procedure is repeated 
by changing signal-to-noise ratio (SNR) of the supernova signals, which is the optimal 
output of the matched filter method. In this algorithm, we did not include the step 
to find the peak positions because SNR of signals is large enough to distinguish the 
difference of the second peak and the third peak and the positions of the peak are 
seen clearly. 

Figure [5] shows the detection probability of the second peak as a function of 
the averaged signal-to-noise ratio over all detectors. From this figure, at 3:00 am, 
December 25th in 2007, the second peak of the signal MCS, RCS with signal-to-noise 
ratio > 6 is detectable with > 90% efficiency. The signal-to- noise ratio of the supernova 
signals at the location of Galactic center for model MCS is : 6.8 for HI, 5.5 for HI, 
13.4 for LI, 9.6 for VI, and for model RCS: 8.9 for HI, 6.4 for H2, 12.7 for LI, 4.9 
for VI. This shows that the detection of the supernova signals at the Galactic center 



Determination of the angular momentum distribution of supernovae 



-4 



-6 



jl/yw^NsWAAAr'U.y'^ ^_..,~~-»s~— 



0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 
time [s] 



Figure 5. The reconstructed waveform of model SCS. The signal-noise-ratio is 
30. 



enables us to infer the dynamics of the supernova core, particularly the rotation law 
and the degree of the differential rotation. Here, signal-to-noise ratio is the output 
obtained from matched filtering and is calculated by using design sensitivity of each 
detector jI^HT]. 

Regarding the signal SCS, the detection probability is reduced as signal-to-noise 
ratio increases in spite of the the fact that the signal-to-noise ratio of model SCS is 3.8 
for HI, 3.0 for H2, 8.2 for LI, 6.8 for VI. When signal-to-noise ratio is small, accidental 
detections occur. This is the reason why detections occur in the region with low signal- 
to-noise ratio. Since accidental detections become smaller and true detections become 
larger as signal-to-noise ratio increases, the detection efficiency becomes smaller. This 
result shows the second peak of the signal SCS cannot detected by the LIGO- VIRGO 
network at this sky location (See figure [5]). 

Figure [6] shows the effect of the regulator on the detection of the second peak. 
A supernova signal with the model MCS was simulated assuming it arrives from the 
direction of SN 1987A. Assuming the event occurs at 0:00 am, December 25th in 2007, 
and other conditions being the same as the previous simulations, we calculated the 
detection efficiency. The red line in the figure is the result obtained by the regularized 
coherent network analysis with regulator gain equal to 1 . The regulator gain is a scalar 
value taking on values in the range [0, 1] and controls the strength of the regulator. 
The blue line is the result obtained by the non-regularized coherent network analysis, 
i.e., the regulator gain being zero. Since the regulator works to reduce variance of 
the extracted signals, it plays a significant role in the astrophysical interpretation of 
detected signals. 
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Figure 6. Detection probability of second peak of the type MCS supernova 
signals at SN1987A as a function of a signal-to-noise ratio. The 'gain' in the plot 
is a scalar value which controls the strength of the regulator, taking values in the 
range [0, 1]. The efficiency obtained by the regularized coherent network analysis 
is significantly superior to that of the non-regularized one. 



5. Conclusion 

From recent siraulations of the dynamics of a supernova core, one finds tliat for 
a cylindrical rotation law with strong differential rotation, the sign of the peak 
with the second largest amplitude is negative. In this paper, wc performed Monte 
Carlo simulations to study the detection efficiency of the sign of the second peak 
using the supernova models MCS, RCS, SCS, and using a coherent network analysis 
method called RIDGE. Using an idealized network of GW detectors, the second peak 
of the signal MCS, RCS with signal-to-noise ratio > 6 is detectable with > 90% 
for the specific source sky location we considered. These results are encouraging 
and indicate that we may be able to obtain information about the dynamics of the 
supernova core, if supernova events occur within our Galaxy. The final answer requires 
refined simulations. We are in progress for this research including a part of the data 
conditioning. We also showed that the regulator imposed on the standard coherent 
network analysis plays an important role in improving the detection efficiency of the 
second peak. 

In figure [1] a low frequency component, which is like a direct current (DC) 
component, can be observed in the waveforms. This component appears in ring-down 
phase. In a more realistic case, this component will affect the detection efficiency 
of the second peak since this component is difficult to see due to seismic noise. 
This component comes from gravitational waves radiated from material undergoing 
convection etc in the core and this component has negative sign [371130] . So we can still 
say the negative second peak indicate strong differential rotation of the core compared 
with signals having a positive second peak. 
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